function showCycloneDesign(cutsize,Q,denP,denF,vis,pressure,Eu,Stk50)
set(gca,'xtick',[]);
set(gca,'ytick',[]);
axis([-1,12,-38,1]);
showValues(cutsize,Q,denP,denF,vis,pressure,Eu,Stk50);
showNumberCyclones(cutsize,Q,denP,denF,vis,pressure,Eu,Stk50);
end

function showValues(cutsize,Q,denP,denF,vis,pressure,Eu,Stk50)
text('FontSize',13,'Position',[0 0],'Interpreter','latex','String','$$\underline{Parameters\hspace{2mm}given\hspace{5mm}}$$');
%size
strSize=strcat('$$Required \hspace{2mm}Cut \hspace{2mm}Size\hspace{2mm}x_{50}=',sprintf('%.4f',cutsize*10^6),'\hspace{1mm}\times\hspace{1mm}10^{-6}\hspace{1mm}m$$');
text('FontSize',13,'Position',[0 -1*1.1],'Interpreter','latex','String',strSize);
%denF
strDenF=strcat('$$Fluid\hspace{2mm}Density\hspace{2mm}\rho_f=',sprintf('%.3f',denF),'\hspace{1mm}kg/m^3$$');
text('FontSize',13,'Position',[0 -2*1.1],'Interpreter','latex','String',strDenF);
%denP
strDenP=strcat('$$Particle\hspace{2mm}Density\hspace{2mm}\rho_p=',sprintf('%.3f',denP),'\hspace{1mm}kg/m^3$$');
text('FontSize',13,'Position',[0 -3*1.1],'Interpreter','latex','String',strDenP);
%vis
strVis=strcat('$$Viscosity\hspace{2mm}\mu=',sprintf('%.6f',vis),'\hspace{1mm}Pa\cdot s$$');
text('FontSize',13,'Position',[0 -4*1.1],'Interpreter','latex','String',strVis);
%gasflowrate
strQ=strcat('$$Fluid\hspace{2mm}Volumetric\hspace{2mm}Flow\hspace{2mm}Rate\hspace{2mm}Q=',sprintf('%.4f\\hspace{1mm}m^3/s$$',Q));
text('FontSize',13,'Position',[0 -5*1.1],'Interpreter','latex','String',strQ);
%allowable pressure
strPressure=strcat('$$Pressure\hspace{2mm}Drop\hspace{2mm}across\hspace{2mm}Cyclone\hspace{2mm}\Delta p=',sprintf('%.4f',pressure),'\hspace{1mm}Pa$$');
text('FontSize',13,'Position',[0 -6*1.1],'Interpreter','latex','String',strPressure);
%Euler
strEu=strcat('$$Euler\hspace{2mm}Number\hspace{2mm}Eu=',sprintf('%.3f',Eu),'$$');
text('FontSize',13,'Position',[0 -7*1.1],'Interpreter','latex','String',strEu);
%Stokes
strStk=strcat('$$Stokes\hspace{2mm}Number\hspace{2mm}Stk_{50}=',sprintf('%.6f',Stk50),'$$');
text('FontSize',13,'Position',[0 -8*1.1],'Interpreter','latex','String',strStk);
end

function showNumberCyclones(cutsize,Q,denP,denF,vis,pressure,Eu,Stk50)
text('FontSize',13,'Position',[0 -10*1.1],'Interpreter','latex','String','$$\underline{Determining\hspace{2mm}Number\hspace{2mm}of\hspace{2mm}Cyclones\hspace{2mm}Required\hspace{2mm}to\hspace{2mm}Achieve\hspace{2mm}Minimum\hspace{2mm}Cut\hspace{2mm}Size\hspace{24mm}}$$');
%v
v=sqrt(2*pressure/Eu/denF);
strV=strcat('$$Characteristic\hspace{2mm}velocity\hspace{2mm}\nu=\sqrt{\frac{2\Delta p}{Eu\rho_f}}\hspace{2mm}=\hspace{2mm}',...
    sprintf('\\sqrt{\\frac{2(%.4f)\\hspace{4mm}}{(%.3f)(%.3f)\\hspace{7mm}}}\\hspace{2mm}=%.5f\\hspace{2mm}m/s$$',pressure,Eu,denF,v));
text('FontSize',13,'Position',[0 -12*1.1],'Interpreter','latex','String',strV);
%D
k=sqrt(4*Q/pi/v);
strD=strcat('$$Diameter\hspace{2mm}D=\sqrt{\frac{4Q}{n\pi\nu\hspace{2mm}}}\hspace{2mm}=\hspace{2mm}',...
    sprintf('\\sqrt{\\frac{4(%.4f)\\hspace{4mm}}{n\\pi(%.5f)\\hspace{7mm}}}\\hspace{2mm}=\\frac{%.5f\\hspace{2mm}}{\\sqrt{n}}$$',Q,v,k));
text('FontSize',13,'Position',[0 -14.5*1.1],'Interpreter','latex','String',strD);
%Stk50
strStk='$$Stk_{50}\hspace{1mm}=\hspace{1mm}\frac{x_{50}^2\rho_p\nu}{18\mu D}$$';
text('FontSize',13,'Position',[0 -16.5*1.1],'Interpreter','latex','String',strStk);
strsub='$$Substituting\hspace{2mm}the\hspace{2mm}values:$$';
text('FontSize',13,'Position',[0 -18.5*1.1],'Interpreter','latex','String',strsub);
strStkCal=sprintf('$$%.6f\\hspace{2mm}=\\hspace{2mm}\\frac{(%.4f\\hspace{1mm}\\times\\hspace{1mm}10^{-6})^2(%.3f)(%.5f)\\hspace{14mm}}{18(%.6f)(%.5f/\\sqrt{n})\\hspace{6mm}}$$',...
    Stk50,cutsize,denP,v,vis,k);
text('FontSize',13,'Position',[0 -20.5*1.1],'Interpreter','latex','String',strStkCal);
%n
D=cutsize^2*denP*v/18/vis/Stk50;
n=(k/D)^2;

strN=sprintf('$$\\rightarrow\\hspace{2mm}n=%.3f\\bf{=%d}$$',n,ceil(n));
text('FontSize',13,'Position',[0 -22.5*1.1],'Interpreter','latex','String',strN);
%% Find diameter
text('FontSize',13,'Position',[0 -24*1.1],'Interpreter','latex','String','$$\underline{Determining\hspace{2mm}Diameter\hspace{2mm}of\hspace{2mm}Cyclone\hspace{12mm}}$$');
%D
D=k/sqrt(ceil(n));
strD=sprintf('$$D=\\frac{%.5f\\hspace{2mm}}{\\sqrt{n}}=\\frac{%.5f\\hspace{2mm}}{\\sqrt{%d}}=\\bf{%.3f\\hspace{1mm}m}$$',k,k,ceil(n),D);
text('FontSize',13,'Position',[0 -26*1.1],'Interpreter','latex','String',strD);
%% Find cut size
x50=sqrt(Stk50*18*vis*D/denP/v);
text('FontSize',13,'Position',[0 -28*1.1],'Interpreter','latex','String','$$\underline{Determining\hspace{2mm}Actual\hspace{2mm}Cut\hspace{2mm}Size\hspace{12mm}}$$');
strx50=strcat('$$x_{50}=\sqrt{\frac{18Stk_{50}\mu D}{\rho_p\nu}}=',...
    sprintf('\\sqrt{\\frac{18(%.6f)(%.6f)(%.3f)\\hspace{14mm}}{%.3f(%.5f)\\hspace{7mm}}}=\\bf{%.4f\\hspace{1mm}\\times\\hspace{1mm}10^{-6}\\hspace{1mm}m}$$',...
    Stk50,vis,D,denP,v,x50*10^6));
text('FontSize',13,'Position',[0 -30*1.1],'Interpreter','latex','String',strx50);
end   